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Abstract 

A pure birth stochastic process with several initial conditions is considered. 
We analyze multiplicity distributions of e + e~ collisions and e-p collisions, usig the 
Modified Negative Binomial Distribution (MNBD) and the Laguerre-type distri- 
bution. Several multiplicity distributions show the same minimum x 2 's values in 
analyses by means of two formulas: In these cases, we find that a parameter N 
contained in the MNBD becomes to be large. Taking large N limit in the MNBD, 
we find that the Laguerre-type distribution can be derived from it. Moreover, from 
the generalized MNBD we can also derive the generalized Glauber-Lachs formula. 

Finally stochastic properties of QCD and multiparticle dynamics are discussed. 

PACS number(s): 13.65.+1, 34.10,+x 



1. Introduction 

The Negative Binomial Distribution (NBD) played the important role in analyses of 
multiplicity distributions in high energy collisions §]. However, to resolve the problem 
concerning a parameter k of the NBD the following scheme adopted from QCD was 
proposed in Refs. |3|, §| . In QCD we can consider a pure gluon bremsstrahlung process, 
g — > g + g. The evolution equation is 

= -XnP(n, t) + X(n - l)P(n - 1, t), (1) 
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where P(n,t), t oc ln[ln(Q 2 / /i 2 ) / ln(Ql/ /i 2 )], and A = iV c (color factor) je (cutoff) are the 
gluon distribution, a proper variable in QCD, and the production rate of gluons, respec- 
tively. Assuming n gluons at t — 0, i.e., <5„ ino , they obtained the following distribution 
with n{t) = n exp(Xt), 

P„.(M) = , ,,, (»(! - (2) 

[n — no)\{no — 1)1 n n 

Using the local parton duality between QCD and multiparticle dynamics Eq.(2) was 
applied to the data of high energy collisions in Refs. || ||. However, Eq.(2) does not 
show satisfactory results, because of the fixed multiplicity at t = 0. This suggests that 
the fluctuation at t = is important. 

Later to improve the defects of Eq. (2), two initial conditions, the Poissonian distri- 
bution and the binomial distribution have been considered: Using Eq. (1) with them, the 
Lagueree-type distribution in Ref. || and the Modified Nagative Binomial Distribution 
(MNBD) in Refs. §-|T§ have been derived. See also Refs.|n|-|g. 



In the present paper we analyze the data of e + e~ collisions including the data of 
e-p collisions in Refs.||14|-||21"| by means of three formulas. Several results show that the 
minimum x 2 values by the MNBD with large N contained in the formula are coincided 
with those by the Laguerre-type distribution. See Table I. Thus we investigate whether 
or not there is any interrelation between the MNBD and the Laguerre-type distribution. 

Contents of this paper are as follows: In the next paragraph and the third one, we 
introduce the MNBD and the Laguerre-type distribution, respectively. In the fourth one, 
we consider the large N limit of the MNBD. In the fifth one, we consider the generalized 
MNBD. Finally our concluding remarks are presented. 

2. MNBD and its application to data of e + e~ and e-p collisions 

We assume that the solution of Eq.(l), i.e., P(n,t), is useful for description of mul- 
tiplicity, i.e., in multiparticle dynamics, using the local parton duality. The MNBD is 
obtained from Eq.(l), using the following initial condition, 



P(n, t = 0)= N C n a n (l - a) N ~ n , (3) 

where N is an integer and Na is the average number of clusters. The generating function 
of the MNBD is 

U(u) = £ P{n)u n = [1 - n{u - l)] N [l - r 2 (u - 1)]~ N , (4) 

n=0 

where r 1 = p — a(p+ 1), and r 2 = p with p = exp(At) — 1. In this case A is the production 
rate of hadrons. 
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The MNBD is given by the following function |], [J, neglecting the second variable, 

"I +7-1-1" 



P(0) 
P(n) 



.l + r 2 J 

1 rl + rnJV 



n! 



X 



i + r j E»c M) 1 

B ns (N + n-s-l)\ 



(N-s)\ 



(N-l)\ 



(5) 



where A = 1 + ri/(l + ri) and 5 = r2/(l + r 2 ). The ri and r 2 are expressed by the 
average multiplicity (n) and the second moment C 2 = (n 2 )/(n) 2 . In this case r\ is real 
and r 2 > 



ri = - C 2 - 1 - - (n) - -, r 2 = n + V- 



AT 



AT 



(6) 



Our analyses of the charged multiplicity distributions in e + e~ and e-p collisions are shown 
in Table I. 

The KNO scaling function of the MNBD has been calculated by the usual procedure, 



z = lim n/ (n) fixed, in a recent paper [|T0| : 

n,(n)—*oo 



ty(z) = lim (n)P(n) 

n,{n)—> 00 



r' 2 



N 



e '2 



The parameters r' x and r' 2 in Eq. (|7|) are given by 



r 2 — r 



/ \ 3 



r' 2 



(7) 



(c 2 -l-^)(n), r 2 = ri 



+ 



AT ' 



because (n) ^> 1. Our results by the KNO scaling functions are given in Table II. 
3. Laguerre-type distribution and its application to data 

Using Eq. (1) with the Poissonian distribution as the initial condition, 

(in \ 

P(n,t= 0) = ,J fc l_jLexp(-(m», (9) 

we obtain the Laguerre-type distribution. (Notice that (to) = 2/(C 2 — 1 + l/(n)) has 
been interpreted as the number of clusters at t — 0.) 



P(0) = exp(-(m)), 
P{n) = (to) 2 exp(— (to)) 



((n) - (m)) n ~ l 



(n) n 



- L 

n 



(i) 



(to) 1 



(n - 1)V (ti)-(to)' 



n=l,2,... (10) 
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It should be stressed that Eq.(lO) was derived in various approaches and applied to 



analyses of data in high energies in Refs. O, 12, 13 



The KNO scaling function of Eq. ( |T0"D is given as 

= ^■e K p[-(m) k (z + l)]xI 1 (2(m) k Vz), (11) 



where (m)k = 2/(6*2 — 1)- Our results by means Eqs. (10) and (11) are also given in 
Table I and II. Therein we can see that several minimum x 2 's values by Eq.(5) with large 
N are coincided with those by Eq.(10). This fact suggests us that the MNBD in large 
N limit probably contains the Laguerre-type distribution. Thus in the next paragraph 
we investigate this possibility. 

4. Large N limit of the MNBD 

We take the large N limit in Eq. @: 

Eq. (12) is interpreted as the distribution of clusters at the initial time t = as Ref. ||. 
The origin is attributed to the Poissonian distribution at t = 0, which can be reduced 
from the binomial distribution as 

N C n a n (l - a) N ~ n la ^ N ^^exp(-iVa). (13) 

Using Not = (n)/(l + r 2 ) = (m), we obtain Eq.(10) as Ref. ||. Using an approximation 
N Cj ~ we obtain Eq.(10) neglecting the contribution with 0(1/N). Thus it it can 

be said that Eq.(5) contains Eq.(10). For the KNO scaling functions, we have Eq.(|TT|) 
taking the large N limit in Eq.(17), and using the following approximation 

1 -M f _1 /™\2, 



lim -L^A—(m)iz) » 7 - r - P / 1 (2(m) fcv /z). (14) 



The present calculation shows that there is an interesting interrelation between solutions 
of Eq (1). Here it is worthwhile to notice that a similar interrelation between the Ising 
model with the free lattice and Eq.(10) was used in Ref. [12]. 

5. Generalized MNBD 

According to calculations of Refs. [2,3], we can consider Eq.(l) with the immigration as 
follows 

dP{n, t)/dt = -XnP(n, t) + A(n - l)P{n - 1, t) - X P{n, t) + X P(n - 1, t), (15) 

where Ao is an immigration rate corresponding to q — > q + g process in QCD[2-3]. (On 
the other hand, in multiparticle dynamics the physical meaning of Ao may be interpreted 
as possible contributions from neutral excited hadrons.) 
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By making use the same initial condition, Eq. (3), we obtain the following generating 
function with k = A /A as 



U(u) = J2 P(n)u n = [1 - fi(« - 1)] N [1 - r 2 (u - 1)] 



-k-N 



(16) 





n=0 


where 






N 



1 "H C2_1 "R 



r 2 



TVf i + (n) 



Notice that Eq. (16) contains the NBD(N = 0) and the MNBD(/c = 0). The probability 
distribution is obtained from H(u) as 



p(o) = n(o) 



(i + f 2 f +fc: 



p(n) - 1 y cfr-ir An 



s=0 



(N-s)\ 



X 



&n _ s (N + k + n-s-l)\ 
(N + k-l)\ 



(17) 



where A = r\j (1 + fi) and B = f 2 / (1 + r 2 ). 

We obtain the following generalized Glauber-Lachs formula in Ref. [1] taking large N 
limit in Eq.(17), 



N 



Eq. (17) P(n) 



(p(n)/ky 



(l+p(n)/k) 



n+k 



exp[ 



(n> 



(1 +p{n)/k) J 



n 1 (i+p<n)/*) ; ' 



(18) 



where 



„ large N{n) f r —7; " 1 (n) (ft) 



, p, andf 2 - fi 
k N 



6. Concluding remarks 

Equations (5) and (6) and the NBD are used in analyses of the data. We have found 
that several minimum \ 2 values are coincided with each other. Referring this fact, we 
have shown that Eq.(10) is derived from Eq.(5) as N is large. Moreover the generalized 
Glauber-Luchs formula is also derived from Eq.(17) as N is large. 

We consider physical correspondences between QCD and multiparticle dynamics in 
Eq.(l). Using Eq.(4), exp(At) = 1 + r 2 , A and t can be estimated by the data. They 
are shown in Fig. 2 and Table III. We obtain an interesting expressions for A and t in 
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multiparticle dynamics. The different interval concerning t is probably attributed to the 
hadronization. More data at \/s=135 GeV are necessary to confirm them. Finally it 
should be noticed that there is a different approach to stochastic equation in QCD in 
Ref. [22]. 
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Table captions 
Table I: 

The minimum values of x 2 's for fitting of discrete distribution to the experimental data 
using the MNBD the NBD and the Laguerre-type distribution, Eq.(lO), respectively. 
The data of (n) and C2 are used. V denotes value of the minimum x 2 ' s - The NBD is 
given by P(n) = T(n + k)/V{n + 1) F(k) (k/(n)) k (1 + k / '(n))-^ n+k \ 

Table II: 

The same as Table I but for fitting of the KNO scaling function obtained from the 
MNBD, Eq. (7), the gamma distribution and the Laguerre-type distribution, Eq.(ll), 
respectively. The gamma distribution is given by ^(z) = k k /T(k) e~ kz z k ~ l . 

Table III: 

Physical correspondences between A and t in Eq.(l) in QCD and mutiparticle dynamics. 

Figure Captions 

Fig.l: 

The result of fitting by NBD and MNBD to the data observed in e + e~ collisions by 
DELPHI CoUab.[14]. 

Fig.2: 

1 + r 2 = exp(Xt) vs y/s. A function exp(\ln(y/s / s )) is chosen among various fittings by 
the MINUIT program: \ 2 — 28.9. We can estimate A and s in multiparticle dynamics. 
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Eq.(l) 


P(n,t) 


A Ao 


interval of t 


QCD 


gluon dis. 


N c N?-l 
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